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Abstract The t] transition form factor is analyzed for the 
first time in both space- and time-like regions at low and in¬ 
termediate energies in a model-independent approach through 
the use of rational approximants. The t] —> experi¬ 

mental data provided by the A2 Collaboration in the very 
low-energy region of the dilelectron invariant mass distribu¬ 
tion allows for the extraction of the most precise up-to-date 
slope and curvature parameters of the form factors as well as 
their values at zero and infinity. The impact of these new re¬ 
sults on the mixing parameters of the rj-rj' system, together 
with the role played by renormalisation dependent effects, 
and on the determination of the VPy couplings from V Pj 
and P ^Vy radiative decays are also discussed. 

Keywords Transition Form Factors, Fade Approximants, 
77 - 77 ' Mixing, Radiative Decays. 

1 Introduction 


complete the time-like region, colliders access to the 
values through the Py annihilation pro¬ 

cesses. The space-like region of the TFFs are accessed in 
colliders by the two-photon-fusion reaction 
where at the moment the measurement of both virtu¬ 
alities is still an experimental challenge. The common prac¬ 
tice is then to extract the TFF when one of the outgoing lep¬ 
tons is tagged and the other is not, that is, the single-tag 
method. The tagged lepton emits a highly off-shell photon 
with transferred momentum q\ = —Q^ and is detected, while 
the other, untagged, is scattered at a small angle with q^ — 0 . 
The form factor extracted from the single-tag experiment is 
thenF>y*y.(-22,0) =Fpfy{Q^). 

At low-momentum transfer, the TFF can be described by 
the expansion 

Fpry{.Q^) = Fpyy{0) {\-bp^+cp%-dp^ + -A , 
\ nip nip Mp J 

( 1 ) 


The pseudoscalar transition form factors (TFFs) describe the 
effect of the strong interaction on the y*y*P vertex, where 
P = n^,ri,rj', rjc ..and is represented by Fpytyt {q\,q 2 ), a 
function of the photon virtualities q\, and q^. From the ex¬ 
perimental point of view, one can study such TFFs from both 
space- and time-like energy regions. The time-like region of 
the TFF can be accessed at meson facilities either through 
the double Dalitz decay processes P which 

give access to both photon virtualities {q\,q\) in the range 
Arn^ < {q\,q 2 ) < ('«/> —2m/)^, or the single Dalitz decay 
processes P —^ l^l^y, which contains a single virtual pho¬ 
ton with transferred momentum in the range Amj <q\< mj,, 
thus simplifying the TFF to Fpy*f{q\,Q) = Fpy*y{q^). To 
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where Fpyy{Q) is the normalization, the low-energy parame¬ 
ters (LEPs) bp, cp and dp are the slope, the curvature and the 
third derivative of the TFF, respectively, and mp is the pseu¬ 
doscalar meson mass. Fpyy{Q) can be obtained either from 
the measured two-photon partial width of the meson P, 


\Fpyy{0)\^ 


6Ati r{P ^ yy) 
(Ana)^ mj, 


( 2 ) 


or, in the case of T] and 77 ', from the prediction of the 
axial anomaly in the chiral limit of QCD. 

In this work we shall focus on the 77 TFF exclusively. Its 
slope parameter has been extensively discussed from both 
theoretical analyses [1, 3-5, 7] and experimental measure¬ 
ments [8-14]. On the theory side. Chiral Perturbation The¬ 
ory (ChPT) predicts bri = 0.51 at = 0.69 GeV^ and for 
sin 9p = —1/3 [1], being /r the renormalisation scale and 9p 







2 


the rj-rj' mixing angle. Other theoretical predictions are [1]: 
bri = 0.53 from Vector Meson Dominance (VMD), brf = 
0.51 from constituent-quark loops (QL), and bri —036 from 
the Brodsky-Lepage (BL) interpolation formula [2]. Recently, 
the slope has been predicted to be b-q — 0.546(9) and bq — 
0.521(2) from a chiral theory with one and two octets of 
vector resonances [3], respectively, bq — 0.60(6)stat(3)sys 
from rational approximants [4], bq — 0 . 62 )^ qq 3 [5] and bq — 
0.57+qqI [6 ] from dispersive analyses, and bq = 0.51 or 
0.54, depending on the data set used as input, from anomaly 
sum rules [7]. With respect to the experimental determina¬ 
tions, the values for the slope are usually obtained after a fit 
to data using a normalized, single-pole term with an associ¬ 
ated mass Ap, i.e. 


FpfriQ^) 


Fpyy{0) 
l + Q^A^ ■ 


(3) 


The results are bq — 0.428(89) from CELLO [9] and bq — 
0.501(38) from CLEO [10], both from space-like data, and 
bq = 0.57(12) from Lepton-G [8], bq = 0.585(51) from 
NA60[11]^ =0.58(11) from A2 [12], and bq =0.68(26) 
from WASA [13], all of them from time-like data. More 
recently, the A2 Collaboration reported bq — 0.59(5) [14], 
the most precise experimental determination up to date. The 
curvature was for the first time reported in Ref. [4] with 
the value Cq = 0.37(10)stat(7)sys- Nothing is yet reported 
about the third derivative of the TEE, dq, although its role 
on hadronic quantities where TFEs are important suggests 
also to look at it (see Ref. [4] for its role on the 7] contri¬ 
bution to the hadronic light-by-light scattering piece to the 
muon (g — 2 )). 

Several attempts to describe the rj TEE are available in 
the literature at present [3, 5, 7, 15-29] but none of them 
tries for a unique description of both space- and time-like 
experimental data, specially at low energies. In Ref. [30], 
it was suggested for the case that a model-independent 
approach to the space-like TEE can be achieved using a se¬ 
quence of rational functions, the Fade Approximants (PAs), 
to fit the data. Later on, in Ref. [4], the same method was 
applied to the 77 and rj' TFEs. More recently, the A2 Col¬ 
laboration reported a new measurement of the 77 e+e 7 
Dalitz decay process with the best statistical accuracy up 
to date [14]. A comparison with different theoretical ap¬ 
proaches was also performed. In particular, the results from 
Ref. [4], based on space-like data, were extrapolated to the 
time-like region and agreed perfectly with their measure¬ 
ment. Triggered by these new A2 results, we explore in the 
present work a combined description of both space- and time¬ 
like regions of the 77 TEE within our method of rational ap¬ 
proximants. This will provide, for the first time, a determina¬ 
tion of the energy dependence of the 77 TEE in both regions 
together with a unified extraction of its LEPs. 


Our approach makes use of PAs as fitting functions to 
all the experimental data. PAs are rational functions P^{Q^) 
(ratio of a polynomial T^iQ^) of order N and a polynomial 
Rm{Q^) of order M) constructed in such a way that they have 
the same Taylor expansion as the function to be approxi¬ 
mated up to order Since PAs are built in 

our case from the unknown low-energy parameters (LEPs) 
of the TEE, once the fit to the experimental data is done, 
the reexpansion of the PAs yields the desired coefficients. 
Being rational functions the PAs are analytic everywhere 
except where the poles are located. Branch cuts cannot in 
principle be described by PAs, however, if the function to 
be approximated is of a certain kind, for instance a Stieltjes 
function, it can be proven mathematically that an infinite or¬ 
der PA is able to reproduce the cut [3 1 ] ^ Another interesting 
issue is the implementation of chiral logarithms of the kind 
log( 2 ^/M^), appearing for instance in chiral expansions at 
next-to-leading order, in the PAs method. These chiral logs 
admit a Taylor expansion which can be seen as an infinite or¬ 
der diagonal PA and is convergent for any value of > 0. 
Therefore, in case the approximated function includes chiral 
logs their effects are incorporated in the PAs to a good extent 
(more precise as the order of the PA increases)^. The advan¬ 
tage of PAs over Taylor expansions is their ability to enlarge 
the domain of convergence. However, to prove the conver¬ 
gence of a given PA sequence is a difficult task and only 
for certain classes of functions this can be done rigorously. 
In practice, the success of PAs in the description of experi¬ 
mental data can only be seen a posteriori in the sense that 
the pattern of convergence can be shown but unfortunately 
not proven mathematically. We refer the interested reader to 
Refs. [32, 33] for details on this technique. 

In this work, we resume our method [4] for fitting the 77 
TEE experimental data after including all the recent avail¬ 
able time-like measurements from 77 -> l^l^y decays (/ = 
e,li). Besides recapitulating the main features of the method 
we will address the following issues: 

- A reevaluation of the systematic errors considered in our 
previous work is demanded by the inclusion of time-like 
data at these low energies. This new set of data being 
more precise than the space-like one, its incorporation 
will allow for an improved systematic error associated 
with each element of a given PA sequence and the in¬ 
crease in order of the sequence itself 

- The better description of the low-energy region of the 
TEE allows for an improved determination of its value 
at zero momentum transfer, which is related to the two- 

* In thi.s case, the infinite number of poles and zeros are seen to be 
located along the branch cut with the first pole located at the beginning 
of the branch point. 

^ For values of (£- /nF ~ 1 the relative error between the chiral log and 
an associated, for instance, second-order PA is of the order of the per 
mille. 
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photon decay width of the Tj. The impact of the recent 
measurement of this width by the KLOE collaboration 
[34] and of older measurements based on Primakoff tech¬ 
niques is commented. 

- The role played by high-energy space-like data in view 
of the fact that in such region only BaBar data is avail¬ 
able. Related to this, the existing puzzle between the 
precise mixing scenario derived from the TFF in con¬ 
trast to the measured time-like cross section by BaBar 
at = 112 GeV^ [35] is discussed. The possibility for 
the Belle Collaboration to measure the time-like J] TFF 
is also mentioned; 

- The extraction of the rj-rj' mixing parameters from the 
TFFs and the two-photon decays after discussing the 
role of the renormalization scale dependence of the sin¬ 
glet decay constant Fq. The new results are much better 
constrained with the inclusion of time-like experimental 
data and turn out to be competitive with standard deter¬ 
minations, such as for instance the analysis of V — Py 
{V = p,a),0) andPy 7 decays [36]. 

- The determination of these VPy coupling constants from 
the former mixing parameters and its comparison with 
current experimental values. The effect of OZI-violating 
parameters and higher-order effects is also discussed. 

The paper is organized as follows. In Sect. 2, a reanaly¬ 
sis of the systematic error related to our method when tak¬ 
ing into account both space- and time-like experimental data 
is performed. In Sect. 3, a brief description of the general 
method for extracting the low-energy parameters of the 77 
TFF using rational approximants is presented and then the 
impact on them of both the Tj —7 77 latest measurement and 
the high-energy space-like data is discussed. In Sect. 4, the 
implications of our new results for the determination of the 
Tj-Tj' mixing parameters, the understanding of the BaBar 
puzzle, and the prediction of the VPy couplings are exam¬ 
ined. Finally, in Sect. 5 the conclusions of the present anal¬ 
ysis are given. 

2 A new systematic error 

In the context of Fade approximants, by systematic error it 
is meant the difference between the function to be approxi¬ 
mated and the highest approximant reached after the fit pro¬ 
cedure. If there is seen convergence, the larger the PA order, 
the smaller the systematic error. Therefore, any finite order 
PA should have a definite systematic error. In this section, 
we discuss how to obtain such an error for a scenario con¬ 
taining both time- and space-like data. 

In order to illustrate the utility of the PA as fitting func¬ 
tions, Ref [30] simulates the real situation of the experimen¬ 
tal data on the space-like region by generating with differ¬ 
ent models a set of pseudodata. Such data were then fitted 


with a Pi{Q^) (single-pole approximants) sequence and the 
FEPs where extracted. This exercise was twofold: first it was 
meant to show the ability of the PA sequence to extract the 
FEPs and, second, also provided a systematic error for the 
extraction of each FEP at each value of L. In Ref. [4], more 
examples were worked out and further discussed and we re¬ 
fer the interested reader to such references. 

Dealing now with a larger set of data, such systematic 
errors should be reanalyzed, specially because the amount 
of time-like data, which covers the lowest energy region — 
and is most important for FEPs extraction— is larger than 
the space-like one. 

Following the strategy presented in Refs. [4, 30], we 
simulate with an holographic model (see Appendix B for 
details on the model together with on the simulation) the sit¬ 
uation of the experimental data from both space- [9, 10, 37] 
and time-like [11, 12, 14] data. The results obtained with the 
holographic model described in Appendix B are collected in 
Table 1 where the relative errors for the first three deriva¬ 
tives for each element on the Pi{Q^) sequence are reported. 
These results are model dependent. Using, instead, the quark 
model considered in Ref. [30], we find faster convergence 
and we reach systematic errors one order of magnitude bet¬ 
ter for the higher PA of the sequence than the holographic 
model. We chose to use the results with the holographic one 
to be on the conservative side. 

The strategy is then to generate pseudodata for both re¬ 
gions trying to emulate the real experimental situation. In 
the space-like region, we evaluated the model at 10 points 
in the region 0.6 < < 2.2 GeV^, 15 points in the re¬ 

gion 2.7 < 2" < 7.5 GeV^, and 9 more points in the region 
9 < 2^ < 34 GeV^. In the time-like region, the model is 
evaluated at 8 points in the region (0.045)^ < 2 ^ < ( 0 . 100)2 
GeV^, 15 points in the region (0.115)2 < 22 < (0.220)2 
GeV^, and 31 more points in the region (0.230)^ < 22 < 
(0.470)2 GeV^. On top of these set of data points we add 
the value of F,jyy(0,0). All these data points have zero er¬ 
ror because we want to obtain a pure systematic error on 
our fitting functions. Notice that the majority of points lie in 
the low-energy region. This simple exercise also prevents us 
against over-fitting problems. The very same study can be 
performed to evaluate the P^{Q^) sequence. The results are, 
however, an order of magnitude better than for the Pi{Q^) 
one (see the comparison in Appendix B). From now on, we 
consider only the systematic errors from the latter to be on 
the conservative side. 

The forthcoming BESIII data on the t] TFF at space-like 
region below 9 GeV2 might demand a reanalysis of our sys¬ 
tematic errors, although we think that including them would 
not really modify our percentages beyond the precision we 
are reporting them in Table 1 . Their data will be, neverthe¬ 
less, crucial to reduce our statistical errors which is by now 
the dominant source. 
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L 

1 

2 

3 

4 

5 

6 

2 

8 

9 

bq 

9.6 

7.0 

4.3 

3.0 

1.8 

1.1 

0.7 

0.4 

0.2 

Cri 

- 

4.0 

4.0 

3.5 

2.7 

2.0 

1.4 

0.8 

0.5 

dn 

- 

- 

22.2 

18.9 

14.6 

11.3 

8.6 

5.9 

4.0 


value Fri^YY = 0.516(18) keV [42] (which is basically dom¬ 
inated by the recent KLOE-2 measurement [34]) in our fits. 


Table 1 Collection of systematic errors (in percentage %) of the first 
three derivatives bj^^c-q and dq of the for a Pi'(Q^) se¬ 

quence fit. 


In passing, we also study what would be the systematic 
error done by a VMD fit to only the time-like data set. From 
the three models considered in Refs. [4, 30], the most con¬ 
servative systematic error found is around 5% (details are 
presented in Appendix B). Notice that when fitting space¬ 
like data with a VMD such error is around 40%. The reason 
of such difference is simple because available time-like data 
is much closer to the origin of energies than the space-like 
one and less sensible to higher-order effects. 


3 Tj transition form factor: a space- and time-like 
description 

To extract the 7] TFF low-energy parameters bq, Cq, and 
dq (slope, curvature, and third derivative respectively) from 
the available data, we start with a Pi{Q^) sequence. How¬ 
ever, according to Ref. [38], the pseudoscalar TFFs behave 
as 1/2^ for which means that, for any value of L, 

one will obtain in principle a good fit only up to a finite value 
of but not for Therefore, it would be desirable 

to incorporate this asymptotic-limit information in the fits to 
Q^FqytY{Q^) by considering also a P^{^) sequence. 

This method, which makes use of experimental data and 
theoretical framework for fitting them, cannot access the 
second Riemann sheet where the resonance poles are sup¬ 
posed to be located [39]. One cannot extract resonance poles 
parameters with such methods, and that poses a word of cau¬ 
tion on the interpretation of fits such as Fq. (3) to relate its 
pole parameters with effective masses. Our method does not 
contain a branch cut and all the analytical structure is built 
to reproduce only the first Riemann sheet. The effective pole 
we obtain should lie outside the range where data are. The 
main advantage of the method of PAs is indeed to provide 
the 2^ dependence of the TFF over the whole space- and 
time-like region up to the first resonance in an easy and sys¬ 
tematic way, without the need of a model for the resonance 
poles appearing in the amplitude [30, 32]. For how to extract 
resonance pole parameters using PA, see Refs. [40, 41]. 

Experimental data from the space-like region is obtained 
from CEFFO, CFEO, and BABar Collaborations [9, 10, 37], 
together with the time-like experimental data from NA60 
and A2 Collaborations [11, 12, 14]. We also include the 


3.1 A remark on experimental systematic errors 

When comparing time-like data results from different col¬ 
laborations it is common to report, together with the exper¬ 
imental data, the result of a fit with a single-pole function 
Eq. (3). Although such data contain only statistical errors, 
systematic errors are incorporated in the result of the fit. 
When using these data in our fits one must incorporate the 
systematic error information into the fitted data^. 

The A2 Collaboration reported in 2011 on the Dalitz 
decay rj —> [12]. Their fit yielded = (1.92 ± 

0.35jra, ± O.Ojyir) GeV^^. Combining both statistical and 
systematic error one obtains A = (1.92zL0.39comb) GeV^^. 
In order to obtain the combined error from a direct fit to the 
published data one can include a new source of error de¬ 
fined in the following way: ^lat + 

for each 2^ datum, with e a percentage. For the A2 2011 
data we find that e = 6 . 8 % will allow us to reproduce, with 
Eq. (3), the combined result A = (1.92 ±0.39 comb) GeV 

At the same time, an analysis of the rj -A 7 Dalitz 

decay by the NA60 Collaboration allowed a determination 
of A with significantly better statistical accuracy. In 2009, 
they reported the value A^^ = (1.95 ± O.lTi^a, ± 0.05 
GeV^^ [11]^, which implies a factor e = 1.9% for out A fi„ai . 

In 2013, the A2 Collaboration reported a new measure¬ 
ment of the same Dalitz decay rj —> e+e^ywith larger statis¬ 
tics with a fitted value A^^ = (1.95 ± 0.15ira; ± 0.10^^^?) 
GeV-2 [14], which leads to e = 4.8%. 

Published space-like data contains both error sources sep¬ 
arately. The exception is the CEFFO Collaboration which 
does not report a systematic error for each bin of data. Only 
a 12% for the two-photon t]-decay channel is reported. Ac¬ 
counting for all the different systematic sources we could 
find in their publication, we ascribe a 12 % of systematic er¬ 
ror for the hadronic rj decay which leads to a 6 % error for 
the global number of events (implying a 12 % of systematic 
error for each bin). We expect that the forthcoming space¬ 
like measurements at BES-III will provide the accurate de¬ 
scription of such energy region and the role of the unknown 
systematic effects in the CEFFO data would not be impor¬ 
tant. 


^We thank Marc Unverzagt for discussions along these lines. 
^Recently, NA60 presented an improved preliminary result, — 
(1.951 ± 0.059i,„, ± 0.042sys,) GeV“^ [43] but the corresponding data 
are not yet published. 
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3.2 Results 

After defining the set of data we will use, we report on our 
results. We start fitting with a sequence. We reach 

L = 1 and we show it in Fig. 1 as a green-dashed line. The 
smaller plot in Fig. lisa zoom into the time-like region. The 
obtained LEPs are collected in Table 2 and shown in Fig. 2 
together with our previous results (empty orange) when only 
space-like data were included in our fits [4]. The stability 
observed for the LEPs with the Pf (2^) sequence is remark¬ 
able, and the impact of the inclusion of time-like data is clear 
since not only allows us to reach higher precision on each PA 
but also to enlarge our PA sequence by 2 elements. The sta¬ 
bility of the result is also clearer and reached earlier, reduces 
our systematic error, and shows the ability of our method to 
extract, for the first time, the LEPs from a combined fit to all 
the available data. The coefficients of the best fitted Pf (2^) 
can be found in Appendix A. 

To reproduce the asymptotic behavior of the TEE, we 
have also considered the P^{Q^) sequence (second row in 
Table 2). The results obtained are in very nice agreement 
with our previous determinations. The best fit is shown as 
black-solid line in Eig. 1. We reach N — 2. Since these ap- 
proximants contain the correct high-energy behavior built- 
in, they can be extrapolated up to infinity (black-dashed line 
in Eig. 1) and then predict the leading 1 /2^ coefficient; 

hm 2%rr(G") = GeV . ( 4 ) 

This prediction, although larger than in our previous work 
[4], still cannot be satisfactorily compared with the BaBar 
time-like measurement at = 112 GeV^, P,jy«y(l 12 GeV^) = 
0.229(30)(8) GeV [35]. The impact of such discrepancy in 
the rj — rj' mixing is discussed in the next section. 

Our combined weighted average results from Table 2, 
taking into account both types of PA sequences, give 

r =0.576(1 l)s,at(4)sys 

C^=0.339(15)stat(5)sys (5) 

[ =0.200(14)„at(18)sys 

where the second error is systematic (around 0.7, 1.5, and 
9% for bp, cp, and dp, respectively, from Table 1). 

Equation (5) can be compared with bri = 0.60(6)stat(3)sys, 
Ct] = 0.37(10)stat(7)sys Using space-like data exclusively [4]. 
As expected, not only statistical results have been improved 
but also systematics, both by an order of magnitude, yield¬ 
ing the most precise slope determination ever. 

Our slope is compared with experimental determinations 
from [8-14] together with theoretical extraction from [1-7] 
in Eig. 3. 

One should notice that all the previous collaborations 
used a VMD model fit to extract the slope. In order to be 


l-loop ChPT 
VMD 
Quark Loop 
Brodsky Legige 

Di.sp. Rel. 
Axial Anom, 
Lepton-G 
CELLO 
CLEO 
NA60 
A2 
WASA 
A2 


6 
T 
S' 

A 
Iff 
ir 
12 ' 
13 

This WorL 


0.2 


0.4 


0.6 


0.8 


1.0 


Fig. 3 Slope determinations for rj TFF from different theoretical (red 
circles) and experimental (blue squares) references discussed in the 
text. Inner error is the statistical one and larger error is the combination 
of statistical and systematic errors. 


consistent when comparing with our results, a systematic er¬ 
ror of about 40% should be added to the experimental deter¬ 
minations based on space-like data [4, 30], and a systematic 
error of about 5% should be added to the experimental de¬ 
terminations based on time-like data (see Appendix B for 
further details). 

When comparing different theoretical extractions of the 
slope of the ij TEE with our result in Eig. 3, we find a pretty 
good agreement with the exception of the results in Ref [3] 
that reported bri = 0.546(9) and bri = 0.521(2) using Res¬ 
onance Chiral Theory with one- or two-octet ansatze. The 
disagreement is between 2 and 5 standard deviations. Ref¬ 
erence [3] uses Resonance Chiral Theory, which is based 
in large-Ac arguments, to extract LEPs. Going from large- 
Nc to Ac = 3 imposes a systematic error [33, 44—46]. Since 
Ref. [3] considered two approximations for fitting the 7] TEE 
(with one and two octets), one could consider the difference 
between them as a way to estimate such error [4, 40, 47]. In 
such a way, the t] TEE slope would read b-q = 0.53(1), at 
2.5 standard deviation from our result. 

Eventually, we want to comment on the effective single¬ 
pole mass determination Ap from Eq. (3). Using bp = mj,/Ap 
and the values in Eq. (5), we obtain Aq = 0.722(7) GeV or 

2 ^ 1.919(39) GeV-2. 

The fits shown in Eig. 1 use the experimental value of the 
two-photon decay width as an experimental datum to be fit¬ 
ted. Such fit could be repeated without including that decay. 
In such a way, we reach again a P/ (2^) and a (2^) as our 
best PA with the advantage now that the value Pqyy(O) is a 
prediction of our fits. We findP,jyy(0)|/,Y =0.250(38) GeV^^ 
for the P/(2^) and P^yy(0)|/;, = 0.248(28) GeV^^ for the 
P^iQ^), which translates into T]]yy|/,Y = 0.43(13) keV and 
Pqyylfi, = 0.42(10) keV, respectively. Comparing with the 
experimental value Fqyylexp = 0.516(18) keV such predic¬ 
tions are at 0.66 and 0.94 standard deviation each. 
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0 . 10 . 20 . 30 . 40 . 


G" [GeV2] 

Fig. 1 r]-TFF best fits. Green-dashed line shows our best Pf (Q“) fit and black line our best fit. Experimental data points in the space-like 

region are from CELLO (red circles) [9], CLEO (purple triangles) [10], and BABar (orange squares) [37] Collaborations. Experimental data points 
in the time-like region are from NA60 (blue stars) [11], A2 2011 (dark-green squares) [12], and A2 2013 (empty-green circles) [14], The inner 
plot shows a zoom into the time-like region. 



?7 TFF 


N 

brj 

Cri 

dr} 



1 

2 

0.575(16) 

0.576(15) 

0.338(22) 

0.340(20) 

0.198(21) 

0.201(19) 

0.6 

0.6 

Final 


0.576(11) 

0.339(15) 

0.200(14) 



Table 2 Low-energy parameters for the T] TFF obtained from the PA fits to experimental data. The first column indicates the type of sequence 
used for the fit and N is its highest order. The last row shows the weighted average result for each LEP. We also present the quality of the fits in 
tenns of x^/DOF (degrees of freedom). Errors are only statistical and symmetrical. 


3.3 The impact of rj —77 measurements 

Our results in Eq. (5) are, by far, the most precise to date. 
Particularly, we believe that the precision achieved for brj 
will be hard to improve even if new data becomes available. 
Nevertheless, the values obtained mildly depend on TJjyy. 
For instance, if we would have used the value measured 
through the Primakoff mechanism omitted in the PDG av¬ 
erage [42] (i.e, = 0.476(62) keV [48]), we would 

find bri = 0.570(13) for the SL-i-TL extraction with a 20% 
larger error, still in nice agreement with our aforementioned 
results. Notice that this result is pretty similar to the one ob¬ 
tained by our fits when the decay into two photons is not 
used in the data set. This fact does not draw a puzzle, every¬ 
thing seems to agree within uncertainties, but may suggest 
to look again for a Primakoff measurement. 

3.4 The role played by high-energy space-like data 

Low-energy parameters are defined at zero momentum trans¬ 
fer. When extracting them from our fits, one would expect 
the low-energy data to dominate. We noticed, however, that 
in order to reach large PA sequences (leading to more pre¬ 


cise extractions), the high-energy data is also important as 
can be seen in Fig. 1. From 5 GeV^ to 35 GeV^ data are 
basically dominated by the BaBar measurement [37] and 
that has a clear implication on the extraction of the asymp¬ 
totic liniQ 2 ^^Q^FriY*y(Q^} value as can be seen in Table 3, 
where the role of data from each collaboration is reported. 
Indeed, a fit exclusively to BaBar data yields similar re¬ 
sults for both slope and asymptotic value than when con¬ 
sidering the full set of space-like data. However, a fit to the 
data from the CELLO [9] Collaboration which range only 
up to 2.23 GeV^ yields much larger asymptotic value (al¬ 
though statistically compatible). Considering only data from 
the CLEO [10] Collaboration, which ranges up to 12.74 GeV^, 
reduces the asymptotic value by about 20 % compared to 
CELLO. The role of BaBar data is then twofold, allowing 
first to reach A = 2 in the P^{Q^) sequence and determin¬ 
ing basically the asymptotic value. 

In view of the puzzle of the TFF between BaBar [49] 
and Belle [50] results, a second experimental measurement 
covering the high-energy region would be very welcome 
here. We find the Belle Collaboration suited for such pur¬ 
pose and we would like to encourage them to go ahead with 
such measurements. 
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Fig. 2 Slope (top-left panel), curvature (top-right panel), and third derivative (bottom panel) predictions for the rj TFF using the Pi(Q^) up to 
L = 1 (blue points). Previous results considering only space-like data from Ref. [4] are also shown (empty-orange squares) as a way to stress the 
role of the time-like data in our fits. Only statistical errors are shown. 
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On the other side, time-like data can also be used to pre¬ 
dict the asymptotic value, even though the range of data is 
much shorter and much closer to — 0. From the three 
sets of time-like data used in our fits, A2-11 [12] and A2- 
13 [14] are based on the rj —3 7 and covers larger range 

of phase space. The NA60 [11] Collaboration, based on the 
J] —■ covers a shorter range but in the higher-energy 

region. The asymptotic values extracted from the difference 
time-like sets of data agree rather well but disagree with the 
results obtained from the space-like data (although overlap¬ 
ping within errors). Whatever the combination of different 
data sets selected, BaBar data always decides on the asymp¬ 
totic value. In passing, we notice that any of the configura¬ 
tions considered so far agrees with the results of the rj TFF 
measurement from BaR4R [35]. 


4 77 transition form factor: applications 

As sated in the introduction, TFF are not also interesting by 
themselves but also for the range of scenarios where they 
play a crucial role. In this section we consider few of such 
applications. 


4.1 Reanalysis of the rj-rj' mixing parameters 

In this subsection we briefly summarize the main elements 
to extract the mixing parameters exclusively from our fits to 
the form factor data. 

As was done in Ref. [4], we analyze rj-rj' mixing us¬ 
ing the quark-flavor basis. In this basis, the t] and rj' decay 
constants are parametrized as 

( \ ^ -FsSm^s\ 

\FgSm<j)^ F,cos(j)s J ’ 

where F^ s are the light-quark and strange pseudoscalar de¬ 
cay constants, respectively, and ^ the related mixing an¬ 
gles. Several phenomenological analyses And ~ 0^, which 
is also supported by large-Ac ChPT calculations where the 
difference between these two angles is seen to be propor¬ 
tional to an OZI-rule violating parameter and hence small 
[36,51]. 

Within this approximation, the asymptotic limits of the 
TFFs take the form 

lim Q^Fr,y*y{Q^) = 2{c^F^ Fc^F^) 

Q2^oo 

= 2 {cqFq cos 0 - CjF; sin 0 ) , 
lim +CsFA) 

= 2 {cqFq sin -F CiFi cos 0 ) , 


(7) 













data range 
(GeV2) 

pHQ^) 


‘{Q^) 

L 

bj] 

N 

bn 

lim FiqY^ 

02^00 

CELLO [9] 

0.62-2.23 

2 

0.48(20) 

1 

0.427(66) 

0.193(30) 

CLEO [10] 

1.73-12.74 

3 

0.73(12) 

1 

0.522(19) 

0.157(5) 

BABar [37] 

4.47-34.38 

4 

0.53(9) 

1 

0.509(14) 

0.162(3) 

CELLO-rCLEO [9, 10] 

0.62-12.74 

3 

0.65(9) 

2 

0.704(87) 

0.25(10) 

SL 

0.62-34.38 

5 

0.58(6) 

2 

0.66(10) 

0.161(24) 

A2-11-I-A2-13 [12, 14] 

-0.212--0.002 

2 

0.475(76) 

1 

0.551(40) 

0.149(11) 

NA60[11] 

-0.221 --0.053 

3 

0.640(77) 

1 

0.582(19) 

0.141(5) 

TL 

-0.221 --0.002 

3 

0.565(87) 

1 

0.576(17) 

0.143(5) 

CELLO [9] -1- TL 

-0.221-2.23 

5 

0.531(39) 

2 

0.533(30) 

0.203(58) 

CELLO-rCLEO [9, 10] -rTL 

-0.221 - 12.74 

6 

0.567(22) 

1 

0.550(13) 

0.152(3) 

A2-ll+A2-n [12, 14] -rSL 

-0.212-34.38 

7 

0.561(35) 

2 

0.569(28) 

0.178(16) 

TL+SL 

-0.221 - 34.38 

7 

0.575(16) 

2 

0.576(15) 

0.177(15) 


Table 3 Role of the different sets of experimental data in determining slope and asymptotic values of the rj TFF. SL refers the the space-like data 
set, i.e., data from CELLO-tCLEO-tBaB^AR [9, 10, 37] Collaborations, and TL refers to the time-like data set, i.e., data from NA60-rA2-11-rA2- 
13 [11, 12, 14] Collaborations. Bold numbers are our final result. No systematic errors included. 


and their normalization at zero (from the chiral anomaly and 
Eq. (2)) 


FfjyyiO) 


fjj'yy(0) 


47r2 j 


1 ! c,q-c.F« \ 


( 8 ) 


with Cq = 5/3 and Cs = V^/3. 

Experimental information provides |E7)7y(0)|exp =0.274 
(5) GeV-' and |E^'yy(0)|exp = 0.344(6) GeV-^ and for the 
asymptotic value of the t] TEE we take the value shown in 
Eq. (4) with symmetrical errors, liinQ 2 ^^Q^Friyty(Q^) = 
0.177(15) GeV. With these values, the mixing parameters 
are predicted to be 


1.07(1), = 1.39(14), 0=39.3(1.2)°, 

(9) 


where 

f \ = f ^8COS08 -Eosin0o\ 

'v F^, f} J [Fs sin 9g Fo cos9o J ’ ^ ’ 

represents the admixture of the 7] and rj' decay constants in 
terms of the octet and singlet one. 

These relations, Eqs. (10), are very useful since, as ob¬ 
served in Ref. [53] and recently discussed in [54], the singlet 
decay constant Fq is renormalisation-scale dependent; 



^ Fo(m)=Eo(mo)(i + ^('^-‘^)) 

= Eo(aio)(1 + 5), 

with po = — ^Nf, Nc the number of colors, Np the num¬ 

ber of active flavors at each scale, and jj. the renormalization- 
scale, with jLlo = 1 GeV a reference point close to the rj' 
mass. 

To include this effect in our results it is convenient to 
work it out in the singlet-octet basis for later on translate 
it into the flavor one using (10). As such, the asymptotic 
behaviour Eqs. (7) shift to 


with Fjc = 92.21(14) MeV [42]. The uncertainties are dom¬ 
inated by the error from the asymptotic value prediction. 

One can translate the mixing parameters obtained in the 
flavor bases into the octet-singlet one by the following 
recipe [52]: 


F/ = 


F^ + 2F} 


F^ = 


^F^+F^ 


.V2Fs 


9g = (j) — arctan(-), 9q = (j) — arctan( 


V2F, 


( 10 ) 


)■ 


lim Q^F^Y>.y{Q^) = 

Q2^oo 

2 (c^(l -F45/5)7^^cos 0 —Ci(l -|-25)E]sin0), 
lim 2%/y*y(2^) = 

2(c^(l -f 45/5)f’gsin0 -l-Cs(l -|-25)ficos0). 

Assuming asymptotic freedom for o:i(/r), the phenomeno¬ 
logical input as{M^) = 0.1185 [42], and the renormalisation 
group equation for o:i(/r), we determine as{lM) = IGeV) = 
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0.48, including up to four loop corrections and threshold ef¬ 
fects for its running^. With such values and Eq. (12) we de¬ 
termine 5 = —0.17. Using (10) to go back to the flavor basis 
we obtain as our final mixing parameters, representing one 
of the main results of this work: 


inputs :U^yy(0),E)j/yy(0),asympTj 
^ F^/F„ = 1.07(2), Fs/F^ = 1.29(16), 0 = 38.3(1.6)°, 

inputs :Friyy{0),F^ /yy(0),asympTj' 

^ = 1.06(1), = 1.63(8), ()) =41.1(0.8)°, 

(14) 


when taking the ri(ri') asymptotic behavior respectively as 
part of the subset of equations to be solved (8,13). We stress 
that corrections from 5 are bigger for the rj' case, as the 
singlet admixture is more important there. Comparing to our 
old results (sec. Ill in [4]), we hnd a better agreement among 
both solutions in (14). As it was explained in Ref. [4], the 
mixing Eqs. are not independent, there is a relation among 
them: 


lim (Frjy*y(Q^)Friyy(0) +F,^Yy(Q^)^ri'yY(0)) — 




(15) 


where in the last step we have used the exact expressions 
(8,13). Numerically, using our 5 = —0.17, one would ob¬ 
tain 0.85^ for the r.h.s. of (15). Our numerical predictions 
for the asymptotic form factors together with their experi¬ 
mental normalization yield 0.89(3)^ for its l.h.s., in nice 
agreement, but this would not be the case without the 5 cor¬ 
rection. 

This result contrast with BaBar determinations, which, 
taking the running from /To up to their scale = 112 GeV^ 
instead of at (resulting in = —0.09), yields 


= 1.10(3) , FJF^ = 0.91(21) , 0=33(4)°, 
= 1.08(2) , F,/F^ = 1.11(23) , 0=37(3)°, 


using again the two subsets of equations in (13). Eq. (15), 
l.h.s., would read 0.98(7)^ where we are neglecting any 
IjQ^ dependence in it. In other words, assuming that ^ = 
112GeV^ plays the role of oo. The r.h.s of (15) using 
results in 0.92^, then compatible with its l.h.s. This com¬ 
parison is somewhat false since we do not assume SbaBar to 
be at oo, otherwise we would have used d = —0.17, and then 
we would have found a contradiction between l.h.s. and r.h.s. 
of (15). This discrepancy is what we call the BaBar puzzle; 
it is depicted in Eig. 4. This figure includes both our mix¬ 
ing results and the BaBar determination at 2^ = 112GeV^ 

^Particular details of the a,, running are irrelevant at the precision we 
are working. 


FqlFn 


1.058 1.063 1.069 1.078 1.794 1.115 1.160 



Fig. 4 Mixing parameters as a function of the 77 TFF asymptotic value 
when the errors coming from the normalization of the TFFs is set 
to zero. For comparison, we show the mixing parameters extracted 
from the measurement of the time-like rj TFF at 17 ^ = 112 GeV^ by 
BaBar Collaboration [35]. This figure exemplifies the puzzle between 
the standard mixing parameters and the BaBar measurement. 


with the corresponding mixing parameters obtained using 
Eqs. (13,8). 

We remark that our results are tied to a]' corrections 
in ( 12 ) and non-negligible systematic effects for the rj' asymp¬ 
totic behavior from our fits. Since this assertion relies on the 
high energy TEE-behavior, where only BaBar data are avail¬ 
able, a second measurement by Belle Collaboration would 
be a very useful crosscheck. 

The mixing parameters obtained with our fits are pre¬ 
cise enough to be competitive with the standard approaches 
with the advantage of using much less input information. 
Eig. 5 compares our results from Eq. (14) (blue squares) with 
well-established phenomenological determinations, the one 
from Eeldmann, Kroll, and Stech (EKS) from Ref. [51, 52], 
and the one from Escribano and Erere (EE) from Ref. [36] 
(updated in Ref. [4]). The agreement among the three ap¬ 
proaches for both Fy and 0 is impressive. Less agreement is 
found for F^. This parameter is more sensible to meson de¬ 
cays where the strange quark plays an important role, such 
as the 0-meson decays. In fact, such decays where included 
in the EE approach but not in the EKS or in the present work. 

Eig. 6 compares different singlet octet determinations. In 
this basis, our results turn out to be 

Fs/F^ = 1.22(11), Fo/F^ = 1.15(6), 

(17) 

08 = -21.3(3.5)°, 00 =-11.3(3.9)°, 

and they can be compared with the EKS and EE as before 
together with the results by Leutwyler (L) from Ref. [53] 
(no errors were given), and the results from Benayoun, Del- 
Buono, and O’Connell (BDO) [55]. Again, the agreement 
between our results and EKS and EE are remarkable, and 
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also in agreement with the results of Leutwyler. Our results 
slightly disagree for 9o and Fq with BDO. The reason is be¬ 
cause in the BDO the OZl violating piece is not set to zero. 
Since such piece mixes with the singlet component of the 
mixing, their Oq and Fq are slightly shifted. 

4.2 A comment on the BaBar high-energy time-like 
measurement 

Our mixing parameters (14) disagree with the ones obtained 
using BaBar time-like result as shown in Fig. 4 which, in¬ 
cluding the results in Figs. 5 and 6 suggests a puzzle between 
the BaBar measurement and the standard phenomenology [4], 
Assuming that the measurement is correct, the difference 
could be explained if the assumption of duality between time- 
and space-like regions at high energy would not be yet valid 
at 112 GeV^ or that the asymptotic limit of Brodsky and 
Lepage is not yet reached at such energies and a systematic 
error is done in assuming duality [35]. 

In Ref. [35], BaBar collaboration studied the process 
—>■ 7 * —> ^7 at the center-of-mass energy = 

10.58 GeV. They measured its cross section and using its 
relation with the TFF, obtained the absolute value of the 
time-like TFF at =-s = -112GeV^, \Q^F^r-riQ^)\ = 
(0.229±0.031) GeV and \Q^Fr^Yy{Q^)\ = (0.251 ±0.021) 
GeV, where statical and systematic uncertainties are added 
in quadrature. 

A kinematic factor Kp with Kp = I — ^ (see [56]) was 
missing in BaBar expressions. This correction leaves BaBar 
published results almost untouched. This small shift together 
with the duality argument [35], results in a prediction of the 
TFF at 2 ^ ^ 1 12 GeV^: 

\Q^Fr,rr{Q%2=+n2Gev^ = (0-231 ±0.031) GeV, 
\Q%Yr(Q%2=+n2G.v^ = (0.254±0.021) GeV. 

One is tempted to include these time-like measurements 
transformed into space-like predictions (18) into our fits, af¬ 
ter assuming that at this high momentum transfer, the duality 
between time- and space-like holds and no extra error should 
be included. For the rj TFF hts, its inclusion will mainly 
modify the asymptotic prediction growing up its value up to 
lii2aQ2^^Q^Friy*y{Q^) = 0.247 GeV, higher than the BaBar 
result, with a good reduced < 1. This, by itself, already 
indicates that at 2 ^ = 112 GeV^ the asymptotic regime is 
not yet reached. Curiously enough, the value of the fit func¬ 
tion at 2 ^ = -112 GeV^ is | 2 %rr(G^)le 2 =+ii 2 Gev 2 = 
0.219 GeV, below (18). Even worse is the prediction of our 
ht function for the time-like counterpart, i.e., \Q^Frfy*y{Q^) 
|g 2 ^_]i 2 Gev 2 — 0-307 GeV. This exercise shows that the as¬ 
sumption of asymptotic regime at 112GeV^ has an error of 
about 15% in our hts, a theoretical error that should be added 


to BaBar results when used in the space-like region. A re¬ 
cent analysis of the pseudoscalar TFF based on perturbative 
corrections [54] concludes that the difference between the 
time- and space-like form factors at \Q^\ = 112 GeV^ can 
be of the order of 5% to 13% for different pseudoscalar dis¬ 
tribution amplitudes, and can be enhanced by Sudakov-type 
corrections (see [57] for details). The Regge model dehned 
in Ref. [4] also suggests a departure from duality of about 
15% to 20% at 12^1 = 112 GeV^. It is, however, difficult 
to calculate that error and hence difficult to ascribe it to the 
BaBar determination. 

Interestingly enough, to check the eventual departure aris¬ 
ing from duality violations, one could artihcially enhance 
BaBar error just to cross-check its order of magnitude. In¬ 
creasing the error in Eq. (18) fromO.031 GeV^ toO.051 GeV^ 
(adding in quadrature a 1.3(7) and rehtting again, we ob¬ 
tain, whit somewhat better x^, that the asymptotic predicted 
value would then be limg 2 ^oo 2 ^^ 7 )y*r( 2 ^) = 0.193 GeV, 
the ht value at 2^ = -112 GeV^ would read 0.187 GeV 
but also our time-like prediction at 2 ^ = 112 GeV^ would 
read 0.199 GeV, essentially satisfying the initial assump¬ 
tion that time- and space-like TEE coincide at 112 GeV^. 
The error we had to artihcially add to reach at that conclu¬ 
sion is around a 20 %, which agrees with our previous state¬ 
ments and also with [54]. Of course, adding this 20% error 
in Eq. (16) solve what we call BaBar puzzle. 

A 15% departure from the asymptotic limit may seem 
too large for that high momentum transfer. Notice [54, 57] 
that due to its nature, TEE are a convolution of a a pertur¬ 
bative hard-scattering amplitude and a gauge-invariant me¬ 
son distribution amplitude (DA) [58] which incorporates the 
nonperturbative dynamics of the QCD bound-state [38]. That 
means that even for large well inside the asymptotic re¬ 
gion, soft scales coming from the Eock decomposition can 
enhance the TEE. These soft corrections depend on the broad¬ 
ness of the DA. At low energies, our fits suggest the typical 
hadronic scale for the rj TEE to be lower than the rj' coun¬ 
terpart. Being the rj' more contaminated by ss content (and 
less from other Eock states), one would expect its hadronic 
scale to be close to the (j) meson mass, around 1 GeV. This is 
in fact what we find, and indicates a narrower DA for the rj', 
dominated by a qq state, explaining at once why the dual¬ 
ity arguments hold better than in the 77 case. This argument 
complements the one discussed in [54] from the perturbative 
study of the TEEs. 

Even larger error should be added to duality arguments 
at lower energies, such as the measurement of the CLEG 
Collaboration of the same cross section but at ^/s = 3.773 
GeV, and forthcoming measurements by the BES-III Col¬ 
laboration at =4.26 GeV. 

Eor all these reasons, we chose not to use these BaBar 
measurements in the time-like region in our fits. 
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Fig. 5 Mixing parameters of the r\ — r\' system in the flavor basis from different references. 
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Fig. 6 Mixing parameters of the r\ — r\' system in the octet-singlet basis from different references. 


4.3 A prediction for the Vf*/couplings 

In this subsection, we extend our analysis to the vector - 
pseudoscalar electromagnetic form factors. In particular, we 
are interested in the couplings of the radiative decays of 
lowest-lying vector mesons into t] or 77 ', i.e., V —?■ 
and of the radiative decays rj' —^Vy, with V = p,( 0 ,(p. 

We follow closely the method presented in Refs. [36, 
59], and make use of the equations in Appendix A in Ref. [36] 
to relate the form factors with the mixing angle and the de¬ 
cay constants in the flavor basis. To account for the ^ — o 
mixing we use (j)v = 3.4°. The form factors, saturated with 
the lowest-lying resonance and then assuming vector meson 
dominance, can be expressed by 

Fvpy{0,0) = ^gvpy, (19) 

niy 

where gypy are the couplings we are interested in, and fy 
are the leptonic decay constants of the vector mesons and 


are determined from the experimental decay rates via 

r{V^e+e-) = '^a^^cl, (20) 

with cy an electric charge factor of the quarks that make 
up the vector, cy = for V = p,(0,<j) re¬ 
spectively. Here By = -|-arctan(l/-\/2). Experimentally 

we find 

fpo = (221.2 ±0.9)MeV, 

/^ = (179.9±3.1)MeV, (21) 

= (239.0±3.8)MeV. 

using r(p ^ e+e-) = 7.04(6) keV, r(a) ^ e+e-) = 0.60(2) 
keV, and r(0 ^ e+g-) = 1.27(4) keV from [42]. 































12 


The couplings in this flavor basis are: 


Sp’ir - 


3mp cos (j) 
47r2/p0 


8 pr\'y - 


3mp sin (j) 
47r2/p0 ’ 


mco ^ COS 0 sin 0 


V 2 fJ ’ 


111(0 / , sin 0 , , . , COS^ 


V2F, i ’ (22) 


/ COS0 sin^ \ 


mp . , sin0 COS0 \ 


where we have assumed = 0^ = 0. Table 4 collects our 
predictions in its second column. Corrections due to ^ 
to these formulae can be found in Appendix A, Eq. (A.5) of 
Ref. [36]. 

The decay widths of P —>^ V 7 and V Pj are 


r{p^vr) 
r{v ^Pr) 


a J f nip — \ ^ 

a 2 (~ \ ^ 


(23) 


The experimental decay widths from [42] allow us to 
extract an experimental value for gvpy, which are collected 
in the last column on Table 4. 

Our predictions compare well with the experimental de¬ 
terminations, see Table 4, specially considering the simplic¬ 
ity of the approach. The differences are always below 2 stan¬ 
dard deviations, excepting the co couplings. Our prediction 
for the ratio of y/f' decays is in that respect remarkable. 

The observed deviations hint towards a somehow over¬ 
simplified approach. Even though our goal is just to show 
the relevance of TEE in other decays, and we do not pre¬ 
tend an exhaustive study of higher-order contributions in our 
scheme, we still want to remark two possible ways to im¬ 
prove our approach. 

On the one hand, the fact that departs from F-^ in 
Eq. (14) may imply a correction through an OZl-violating 
parameter Ai that appears at next-to-leading order in the La- 
grangian of ;i^PT Large-Ac used to define the mixing equa¬ 
tions, = P;r(l + 2 V 1 / 3 ) [52, 60] which in turns imply 7 ^ 
0s, since (j)q — (ps Ai/3. With the result in Eq. (14), we es¬ 
timate Ai ^ 0 . 2 , in agreement with the naive 1 /Nc counting 
(i.e, Ai ~ 1 /Nc ^ 0.3), and then 0^ — 0s ~ 3.8°. 


The ratio Rj/^/ provides direct information on the angle 
0 ^ since 


Py/y = tan2(0^) 





(24) 


with the J meson mass. The experimental Py/ip ra¬ 
tio defined in last the row in Tab. 4, results in 0^ = (38.1 ± 
0.6)°, which implies 0s = (38.1 -f3.8± 1.6)° = (41.9± 1.6)° 
with the error coming from our determination of 0 in Eq. (14) 
Even though both angles are distinguishable, their impact on 
the gvpy is a shift of the form gypy —>■ gvpy /(cos 0 ^ cos 0 s -f 
sin0^sin0s) [36]. Eor the 0,^ = 0s limit, such shift is exactly 
1. Using the 3.8° difference, such shift translates into 0.998, 
a 2 per mil effect, negligible. Our assumption 0^ = 0s = 0 is 
supported phenomenologically. 

On the other hand, as discussed in detail in Ref. [61], in 
the flavor singlet channel one has to allow for another OZI- 
rule violating correction, which essentially corresponds to 
replacing Pq —> Po/(l +A 3 ). This shifts both the P —> 77 
decays and the formulae for gvpy predictions [60]. The pa¬ 
rameter A 3 is, however, still unknown, although expected to 
be ~ 1 /Nc ^ 0.3. We can make use of Eq. (15) to estimate it. 
The shift on Pq can be translated into a shift in P^ s recalling 
that both are related to Fj( following Eqs. (17,9) respectively, 
and find P^ s —> ^ 97 /(1 =*= 21 . 3 ) as well. Going then to Eq. ( 8 ), 

-^77(')77(0) -f77(')77(0)(l ~2V3). 

Then, Eq. (15) transforms into 


lim Q (Friy*y(Q )Friyy(0) + Fp/y*y(Q )Fplyy(0)^ 

Q^^OO 

x(1-^3)=(i + Y)^- 


(25) 


which after expanding and reorganizing in such a way that 
in the l.h.s. remains only experimental quantities, results in: 


lim (Fpy*y(Q^)Fpyy(0) + Fpfy,y(Q^)Fp,yy(0)) 


= (^l-fY +^3 + Y^ 3 ) 


3 

27r2 ■ 


(26) 


We recall that l.h.s., experimentally, reads 0.89^, and 
d = —0.17. With (26) we find A 3 = 0.05, smaller than ex¬ 
pected and with positive sign. 

The VP 7 couplings are also shifted by A 3 . The expres¬ 
sions can be found in Eq. (42) in Ref. [60] which, after ex¬ 
panding, can be expressed as a shift on the couplings in our 
Eq. (22): gVrjy —>■ (gVriy+ \gVTiy\-^3/2) and gvTj'y (^Vt ]'7 + 
\gvr}'y\-^3)’ always increasing the coupling. Eor some of them, 
the A 3 correction goes on the right direction (the p case), but 
for others it is not conclusive (the 0 case where for rj goes 
well and for Tj' wrong). The result of the shift is, then, am¬ 
biguous. 
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Discarding OZI-violating effects, Fade approximants can 
then be the avenue to follow since the vector mass that should 
be used in Eq.(19) it should not correspond to a physical ob¬ 
servable, but an effective scale provided by the pole of a PA 
assuming the philosophy of the present work. For the J] TFF, 
the from (3) is smaller than the vector meson dominance 
mediator. If the same would happen with the p, (O form fac¬ 
tors, one would expect, then, different gvpy couplings. Since 
this study is beyond the scope of the present analysis, we 
postponed for future work. A naive estimate of these effects 
could be accounted for within the half-width-rule [46], i.e., 
instead of using niy in (19), use my ±1^/2, with F the full 
width of the vector. This provides a way to asses the error of 
neglecting the width of the resonance in using my. For ex¬ 
ample, for the p case, within the half-width-rule, the errors 
of the gppy would be enlarged by a factor 3, well compatible 
with the experimental determinations. 

Further studies along these lines are postponed for future 
work. 


Prediction Experiment 


gprjr 

1.50(4) 

1.58(5) 

gpri'r 

1.18(5) 

1.32(3) 

ScoriY 

0.57(2) 

0.45(2) 

ScoTj'y 

0.55(2) 

0.43(2) 

gprir 

-0.83(11) 

-0.69(1) 


0.98(14) 

0.72(1) 

„ _ r(J/V^r,'y) 

4.74(55) 

4.67(20) 


Table 4 Summary of VPy couplings. Experimental determinations are 
from Ref. [42]. 


5 Conclusions 

In the present work, the rj transition form factor has been an¬ 
alyzed for the first time in both space- and time-like regions 
at low and intermediate energies making use of a model- 
independent approach based on the use of rational approx¬ 
imants of Fade type. The model independence of our ap¬ 
proach is achieved trough a detailed and conservative evalu¬ 
ation of the systematic error associated to it. The new set of 
experimental data on the J] -4 7 reaction provided by 

the A2 Collaboration in the very low-energy part of the time¬ 
like region allows for a much better determination of the 
slope and curvature parameters of the form factor, as com¬ 
pared to the predictions obtained in our previous work only 
using space-like data, which constitute the most precise val¬ 
ues up-to-date of these low-energy parameters. Our method 
is also able to predict for the first time the third derivative of 
the form factor. In addition, the new analysis has served to 


further constrain its values at zero momentum transfer and 
infinity. We have seen that our results, in particular for the 
case of the slope parameter, are quite insensitive to the val¬ 
ues used in the fits for the two-photon decay width of the 
77 , thus showing that the collection of space- and time-like 
experimental data is more than enough to fix a value for 
the normalization of the form factor compatible with cur¬ 
rent measurements. We have also seen that the role played 
by the high-energy space-like data is crucial to get accurate 
predictions for the low-energy parameters of the form factor 
and its asymptotic value. As a consequence of these new re¬ 
sults, we have fully reanalyzed the rj-rj' mixing parameters 
this time also considering renormalisation-scale dependent 
effects of the singlet decay constant Fq. The new values ob¬ 
tained are already competitive with standard results having 
the advantage of requiring much less input information. Re¬ 
lated to this, we have also obtained predictions for the VPy 
couplings which are in the ballpark of present-day determi¬ 
nations. 

In summary, the method of Fade approximants has been 
shown to be very powerful for fixing the low-energy proper¬ 
ties of the J] transition form factor making their predictions 
more accurate and well-established. This fact opens the door 
to a more exhaustive analysis of the single Dalitz decay pro¬ 
cesses P —> l^Fy, with P — 77 *^, 77 , 77 ' and I = e,p, the dou¬ 
ble Dalitz ones P —^ FFI^F (in all possible kinematically 
allowed configurations) [62], and the rare lepton-pair decays 
P^ l+F —see the e+e^ application in Ref. [63], 
which are usually discussed only in terms of monopole ap¬ 
proximations. Indeed, when this work was being concluded 
the BFSIII Collaboration reported a first observation of the 
77 ' —> 7 process measuring the branching ratio and ex¬ 

tracting the 77 ' transition form factor [64]. This new mea¬ 
surement may put our approach with its back to the wall. 
However, a very preliminary analysis of this recent data in 
comparison with our prediction for this form factor in the 
time-like region exhibits a nice agreement but reveals the ne¬ 
cessity of going beyond the vector meson dominance model 
used in the experimental analysis [65]. 

Appendix A: Best Fade approximant fit 
parameterisation 

In this appendix we provide the parameterizations of our 
best P\'{^) fit for the ^Friyy{Q^). Defining P^iQ^) as; 

Ri(e2) i+ne2 ’ 

the corresponding fitted coefficients® for the Q^FrfY*y{Q^) 
are collected in Table 5. 

^For full precision of the coefficients together with the correlation ma¬ 
trix, contact the corresponding authors. 
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p-TEE 

h 

0.27349 

tl 

1.1771 • 10^^ 

h 

-1.1048-10-3 

U 

2.8861 • 10-3 

H 

2.2974-10-® 


-1.5096-10-’ 

tl 

2.3655-10-'* 

n 

1.9584 


Table 5 Fitted coefficients for our best P/ {Q^) for the Q^Friyj(Q^). 


With the coefficients in Table 5, one can extract the slope 
of the TFF by expanding (A. 1) and normalizing the result as 
inEq. ( 1 ): 

= (ft-nAt = 0.5749 (A.2) 

with niri = 0.547853 GeV, to be compared with the third 
column in Table 2. 


Appendix B; Convergence of the Fade approximant 
sequence 


To test how fast the convergence of our PA sequence is we 
analyze here a simple holographic conhning model presented 
in [18] (and also explored in Ref [17]), based on light-front 
holographic QCD where the correct small behavior (in 
order to simulate conhnement) is introduced using the dressed 
current (see [18] for details)^. 

In this context, the TFF is dehned as (assuming, for sim¬ 
plicity, Tj ~ rjs) 


FriY‘y{Q^) 


A dx 

VA;r2/jt 70 {\+xY 


(B.3) 


where Pqq is the probability of finding the qq component in 
the 77 light-front wave function, and we impose Pqq — 0.5 
for numerics. 

Once the model is dehned, by generating a set of pseu¬ 
dodata as it is done is Sec. 2 we can test how fast the PA 
sequence converge to Q^PrjfyiQ^)- We ht these pseudodata 
with both Pi{QY {QY sequences going up to L = 8 

and N = 4 and we show the convergence pattern for the value 
of the 7] TFF at the origin, Frjyy(O), and its three hrst deriva¬ 
tives (briyCrj, and c/rj) in Fig. 7 (the black square shows the 
corresponding value from Eq. (B.3)). This exercises com¬ 
plements the one studied in the Appendix of Ref. [4]. The 
hrst element of the sequence, the P/ (Q^) contains only two 
parameters and then only Frjyy(O) and brj can be directly 
extracted from the ht. By reexpanding it one could predict 
all the other derivatives but we only show ht outputs. The 

^We do not consider higher-twist components here to keep the model 
easy to use. 


same applies for the second element of the sequence, the 
Pi(Q^) its third derivative. The reader should notice 
the hierarchy pattern of convergence. While Frjyy(O) is ap¬ 
proach very fast and with the P^i^) the error is about a 
2 %, the derivatives are worse predict, being the third one 
the worst. Comparing each prediction with its counterpart 
from the model gives and idea of the systematic error done 
by such procedure (see Ref. [4] for details about the sys¬ 
tematic error from our method). The P^i^) sequence con¬ 
verges much faster, reaching the 10 ^^ with its second ele¬ 
ment. The inconvenience is its growing two-by-two inputs 
for each new element, making the sequence more difficult to 
be performed in fits to real data. These results are model de¬ 
pendent and other models may yield different systematic er¬ 
rors. The trend shown in Eig. 7 is, however, general although 
models such us the logarithmic one studied in Ref. [30] con¬ 
verge faster (reaching smaller errors) even though the sys¬ 
tematic error for the hrst elements on the larger. 
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